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Abstract 

Sphingobium chlorophenolicum Strain L-1 can mineralize the toxic pesticide pentachlorophenol (PCP). We have sequenced 
the genome of 5. chlorophenolicum Strain L-1 . The genome consists of a primary chromosome that encodes most of the 
genes for core processes, a secondary chromosome that encodes primarily genes that appear to be involved in environmental 
adaptation, and a small plasmid. The genes responsible for degradation of PCP are found on chromosome 2. We have 
compared the genomes of 5. chlorophenolicum Strain L-1 and Sphingobium japonicum, a closely related Sphingomonad that 
degrades lindane. Our analysis suggests that the genes encoding the first three enzymes in the PCP degradation pathway 
were acquired via two different horizontal gene transfer events, and the genes encoding the final two enzymes in the 
pathway were acquired from the most recent common ancestor of these two bacteria. 

Key words: horizontal gene transfer, biodegradation, enzyme evolution, pentachlorophenol hydroxylase, tetrachlorohy- 
droquinone dehalogenase, tetrachlorobenzoquinone reductase. 



Introduction 

In the early 20th century consumers embraced the concept 
of "better living through chemistry." Massive quantities of 
chlorinated anthropogenic chemicals such as dichlorodiphe- 
nyltrichloroethane(DDT), lindane, pentachlorophenol (PCP), 
atrazine, tetrachloroethylene, and polychlorinated biphen- 
yls (PCBs) were introduced into the environment with little 
concern over their ultimate fate. Many of these compounds 
persist in the environment and have adverse effects on eco- 
systems. Use of some of these compounds (e.g., DDT and 
PCP) is now restricted or prohibited in the United States, 
although some are still used in developing countries. 

The introduction of anthropogenic chemicals into the 
environment imposes selective pressures that can foster 



evolution of novel pathways that allow microbes to access 
new sources of carbon, nitrogen, or phosphorus and/or to 
detoxify toxic compounds. However, degradation of anthro- 
pogenic chemicals is often inefficient because microbes have 
not yet evolved enzymes that efficiently catalyze the steps 
needed to convert anthropogenic compounds into metabo- 
lites in central carbon metabolism (Copley 2009). 

The pesticide PCP is listed as a Priority Pollutant by the US 
EPA due to its toxicity and persistence in the environment. 
A number of microbes capable of mineralizing PCP have been 
isolated from contaminated sites (Saber and Crawford 1985; 
Apajalahti and Salkinoja-Salonen 1987; Schenk et al. 1990; 
Uotila et al. 1 991 ; Radehaus 1 992; Karn et al. 201 0). The best 
studied of these is Sphingobium chlorophenolicum (Saber 
and Crawford 1985; Steiert and Crawford 1986; Crawford 
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Fig. 1. — PCP degradation pathway. PcpB, pentachlorophenol hydroxylase; PcpD, tetrachlorobenzoquinone reductase; PcpC, tetrachlorohydro- 
quinone dehalogenase; PcpA, 2,6-dichlorohydroquinone dehalogenase; PcpE, maleylacetate reductase; HGT, horizontal gene transfer. 



and Ederer 1999; Takeuchi et al. 2001). Sphingobium chlor- 
ophenolicum has patched together a poorly functioning 
pathway (see fig. 1) that allows mineralization of PCP, 
although degradation is slow and the bacterium is unable 
to grow at high concentrations of PCP (Dai and Copley 
2004). The first three enzymes in the pathway are unusually 
ineffective. PCP hydroxylase (PcpB) has a very low k cat of 
0.02 s~ 1 and shows substantial uncoupling; the enzyme fre- 
quently fails to hydroxylate the substrate and instead releases 
H 2 0 2 as the C4a-hydroperoxyflavin intermediate decomposes 
(Hlouchova K, Rudolph J, Pietari JMH, Behlen LS, Copley SD, 
unpublished data). Tetrachlorobenzoquinone (TCBQ) reduc- 
tase (PcpD) is also rather inefficient, with a turnover rate of 
only 0.7 s" 1 at 50 \xhA TCBQ (Dai et al. 2003). Tetrachlorohy- 
droquinone (TCHQ) dehalogenase (PcpC) is subject to pro- 
found substrate inhibition; binding of its aromatic substrates 
to the active site before completion of the catalytic cycle inter- 
feres with the final step required to regenerate the free enzyme 
(Warner and Copley 2007a). These findings are consistent with 
the hypothesis that the enzymes in the PCP degradation path- 
way have not yet evolved to the level of function typical of 
most metabolic enzymes. Thus, 5. chlorophenolicum provides 
a window on the process of assembly of a new pathway at 
a relatively early stage in its evolution. 

We recently sequenced the genome of 5. chloropheno- 
licum L-1. Analysis of this sequence and comparison with 
the sequence of the closely related Sphingobium japonicum, 
which degrades lindane (Nagata et al. 2007), provides 
insights into the origins of the PCP degradation enzymes. 
The rDNA genes of S. chlorophenolicum and S. japonicum 
share 97% identity. The phylogenetic distance between 
5. chlorophenolicum and 5. japonicum is ideal to facilitate 



identification of genes that were present in the most recent 
common ancestor of these two Sphingomonads, as proteins 
involved in core processes show >80% pairwise identities at 
the amino acid level. 

Our analysis suggests that the first three enzymes in the 
pathway were acquired by 5. chlorophenolicum by horizon- 
tal gene transfer (HGT) after it diverged from 5. japonicum. 
In contrast, the last two enzymes in the pathway were pres- 
ent in the most recent common ancestor of S. chloropheno- 
licum and S. japonicum. None of the genes encoding the 
PCP degradation enzymes arose by recent duplication 
and divergence of genes within 5. chlorophenolicum. The 
genes occur in two disparate parts of the genome and have 
not yet been integrated into a compact and consistently 
regulated operon. 

Materials and Methods 

Isolation of Genomic DNA from 5. chlorophenolicum L-1 

A sample of 5. chlorophenolicum L-1 was originally depos- 
ited to the American type Culture Collection and designated 
ATCC 39723. The initial characterization of the enzymes in 
the PCP degradation pathway was done with the ATCC 
39723 strain (Orser and Lange 1994; Copley 2000; Dai 
et al. 2003; Warner and Copley 2007b). However, the ATCC 
39723 strain lost the ability to degrade PCP, so a second 
sample of 5. chlorophenolicum L-1 was deposited and 
designated ATCC 53874. Genomic DNA was isolated from 
this strain. 

A 5 ml culture of ATCC medium 1687 {Flavobacterium 
medium) containing 5 |iM PCP was inoculated with a colony 
from a fresh agar plate of half-strength tryptone soy broth 
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containing 50 jiM PCP and incubated with shaking at 30 °C 
overnight. One milliliter aliquots of cells were harvested 
by centrifugation at 16,100 x g for 1 min. Genomic DNA 
was isolated using the Aquapure Genomic DNA kit (Bio- 
Rad, Hercules, CA). Cell pellets from 1 ml of culture were 
resuspended in 300 jllI Genomic DNA Lysis Solution. The so- 
lution was incubated at 80 °C for 5 min to lyse the cells. 
A solution of RNase A (1 . 5 jil, 4 mg/ml) was added, and tube 
was inverted 25 times before incubation at 37 °C for 45 min. 
The temperature of the sample was adjusted to room tem- 
perature and 100 jlxI of Protein Precipitation Solution was 
added. The sample was vortexed for 20 s and then subjected 
to centrifugation at 1 3,000 x g for 3 min to remove cellular 
debris. The supernatant containing DNA was transferred to 
a clean 1.5 ml tube, incubated on ice for 5 min, and then 
subjected to centrifugation again. The supernatant was 
then transferred to a clean 1.5-ml tube containing 300 jllI 
100% isopropanol to precipitate the DNA. The tube was 
inverted 50 times. The DNA was pelleted by centrifugation 
at 13,000 x g for 1 min. The supernatant was poured off, 
and the pellet was washed with 300 uJ 70% ethanol. The 
ethanol was poured off and residual ethanol removed by 
draining the tube onto filter paper for 1 0-1 5 min. The pellet 
was dissolved in 100 uJ 10 mM Tris-HCI, pH 7.5, at65°C for 
5 min and then at room temperature overnight. The quality 
of the DNA was checked on a 0.8% agarose 1X TAE 
gel alongside a 1 kb DNA Extension Ladder (Invitrogen, 
Carlsbad, CA), and the concentration was measured on 
a NanoDrop spectrophotometer. 

Genome Sequencing 

The draft genome of 5. chlorophenolicum L-1 was gener- 
ated at the DOE Joint genome Institute (JGI) using a combi- 
nation of lllumina (Bennett 2004) and 454 technologies 
(Margulies et al. 2005). Three libraries were constructed 
and sequenced: an lllumina GAii shotgun library that gen- 
erated 40,283, 1 93 reads totaling 3,061 Mb; a 454 Titanium 
standard library that generated 302,660 reads and a paired 
end 454 library with an average insert size of 1 1 .26 ± 2.81 
kb that generated 1 74,361 reads. In total, 1 76.8 Mb of 454 
data were obtained. General aspects of library construction 
and sequencing performed at the JGI can be found at http:// 
www.jgi.doe.gov/. The initial draft assembly contained 79 
contigs in one scaffold. The 454 Titanium standard data 
and the 454 paired end data were assembled with Newbler, 
version 2.3. The Newbler consensus sequences were com- 
putationally shredded into 2 kb overlapping fake reads 
(shreds), lllumina sequencing data was assembled with VEL- 
VET, version 0.7.63 (Zerbino 2008), and the consensus 
sequences were computationally shredded into 1.5 kb 
shreds. The 454 Newbler consensus shreds, the lllumina 
VELVET consensus shreds, and the read pairs in the 454 
paired end library were integrated using parallel phrap, 



version SPS — 4.24 (High Performance Software, LLC). The 
software Consed (Ewing and Green 1998; Ewing et al. 
1 998; Gordon et al. 1 998) was used in the finishing process, 
lllumina data were used to correct potential base errors and 
increase consensus quality using the software Polisher devel- 
oped at JGI (Lapidus A, unpublished data). Possible misas- 
semblies were corrected using gapResolution (Han C, 
unpublished data) or Dupfinisher (Han 2006) or by sequenc- 
ing cloned bridging PCR fragments. Gaps between contigs 
were closed by editing in Consed, by PCR, and by Bubble 
PCR (Cheng J-F, unpublished data) primer walks. A total 
of 87 additional finishing reactions (either sequencing of 
bridging PCR fragments or primer walking) were necessary 
to close gaps and to raise the quality of the finished 
sequence. The total size of the genome is 4,573,221 bp. 
The final assembly is based on 1 76.8 Mb of 454 draft data 
that provides an average 36x coverage of the genome and 
3,553 Mb of lllumina draft data that provides an average 
790 x coverage of the genome. 

Genes were identified using Prodigal (Hyatt etal. 2010) as 
part of the Oak Ridge National Laboratory genome annota- 
tion pipeline, followed by a round of manual curation using 
the JGI GenePRIMP pipeline (Pati et al. 201 0). The predicted 
coding sequences were translated and used to search the 
National Center for Biotechnology Information nonredun- 
dant database, UniProt, TIGRFam, Pfam, PRIAM, KEGG, 
COG, and InterPro databases. These data sources were com- 
bined to assert a product description for each predicted pro- 
tein. Noncoding genes and miscellaneous features were 
predicted using tRNAscan-SE (Lowe and Eddy 1 997), RNAm- 
mer (Lagesen et al. 2007), Rfam (Griffiths-Jones et al. 2003), 
TMHMM (Krogh et al. 2001), and signalP (Bendtsen et al. 
2004). 

Results and Discussion 

The 5. chlorophenolicum Genome 

The 5. chlorophenolicum genome consists of two chromo- 
somes and a plasmid (see fig. 2). A summary of the charac- 
teristics of each replicon is provided in table 1. Like many 
bacteria such as Burkholderia pseudomallei (Holden et al. 
2004) and Vibrio cholerae (Heidelberg et al. 2000) that con- 
tain multiple replicons, 5. chlorophenolicum has a dominant 
chromosome that carries most of the genes for core pro- 
cesses. We generated a list of genes involved in replication, 
transcription, translation, cell division, peptidoglycan bio- 
synthesis, and core metabolic processes (including glycoly- 
sis, the pentose phosphate pathway, the TCA cycle, electron 
transport, and biosynthesis of amino acids, nucleotides, and 
cofactors) (see supplementary table 1, Supplementary Ma- 
terial Online). Figure 3 shows the distribution of genes be- 
tween the two chromosomes. If genes encoding core 
functions were randomly distributed between the chromo- 
somes, we would expect 354 (72%) to be on chromosome 1 
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Fig. 2. — Circle diagram of the replicons of Sphingobium chlorophenolicum. The second circle in each replicon indicates the locations of PCP 
degradation genes in black, phage genes in blue, and transposon-related genes in red. The third circle shows locations of genes that have a close 
homolog in Sphingobium japonicum (>80% identity over >90% of the 5. chlorophenolicum sequence). The fourth circle shows GC content. Red 
indicates sequences with GC content <64% and blue indicates sequences with GC content >64%. Light gray lines are placed at intervals of one 
standard deviation from the mean of 63.8%. (One standard deviation = 0.027%; GC content was calculated by averaging over a 10-kb window and 
sliding that window in 1-kb increments.) Diagram was generated using Circos (Krzywinski et al. 2009). 



and 138 (28%) to be on chromosome 2. In fact, 442 core 
genes (90%) are on chromosome 1 and only 50 (10%) are 
on chromosome 2. This deviation from the expected distri- 
bution is highly significant by chi-square analysis (P value = 
2.2 x 10" 16 ). An odd collection of essential enzymes, includ- 
ing a few genes for central carbon metabolism, two subunits 
of RNA polymerase, some components of the electron trans- 
fer chain, and a few genes for cofactor biosynthesis are 
found on chromosome 2 rather than on chromosome 1 . 
However, most of the genes on chromosome 2 appear to 
be involved in environmental adaptation. Chromosome 2 
carries a number of genes predicted to encode transporters 
for sugars and metal ions and enzymes involved in degrada- 
tion of various sugars, short-chain fatty acids, and aromatic 
compounds. All the PCP degradation genes are found 
on chromosome 2. A secondary chromosome with a low 



density of essential genes may serve as a convenient stor- 
age depot for genes acquired by HGT, as integration of 
newly acquired DNA is not likely to disrupt an essential func- 
tion. Chromosome 2 also contains two genes encoding 
proteins related to ParB, which is involved in plasmid parti- 
tioning. These features of chromosome 2 are consistent 
with the proposal that bacterial secondary chromosomes 
have arisen by intragenomic transfer of essential genes, 
including rRNA genes, to a plasmid (Slater et al. 2009). 

HGT is rampant among microbes and is known to play 
a major role in acquisition of resistance to or degradation 
of toxic compounds, including antibiotics. We examined 
the 5. chlorophenolicum genome for features indicative 
of mobile genetic elements. As mentioned above, 5. chlor- 
ophenolicum contains one plasmid. The genome appears to 
contain one integrated prophage; a cluster of phage-related 



Table 1 

Genome Statistics 





Length (bp) 


GC Content 


ORFs 


rRNA Operons 


Chromosome 1 (NC_015593) 


3080818 


0.64 


2,940 


1 


Chromosome 2 (NC_015594) 


1368670 


0.64 


1,159 


2 


Plasmid (pSPHCHOI) (NC_015595) 


123733 


0.65 


125 


0 
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Fig. 3. — Distribution of genes between chromosomes 1 and 2 of 
Sphingobium chlorophenolicum. Black bars, fraction of all genes 
(excluding those on pSPHCHOI); gray bars; fraction of core genes. 

genes (including a major capsid protein, major tail protein, 
phage portal protein, pro-head peptidase, and some con- 
served phage proteins of unknown function) is found on chro- 
mosome 1 (genes Sc_00004260-00004380). Curiously, seven 
isolated genes annotated as "phage integrase family" 
genes are found on both chromosomes 1 and 2 (Chr1: 
Sc_00026840, Sc_00029370, Sc_00022480, Sc_00028520, 
Sc_00012030; Chr 2: Sc_00031410, Sc_00030440). These 
genes have anomalously low GC content (0.48-0.58) 
and are not closely related to each other. Only one 
(Sc_00026840) is found in the vicinity of a prophage gene 
(a CP4-57 regulatory protein). Two (Sc_00017260 and 
Sc_0003440) are adjacent to transposase genes. The genome 
also carries 27 sequences annotated as "transposase," 
"transposase/integrase core domain," or "transposase and 
inactivated derivatives" (16 on chromosome 1, 9 on chro- 
mosome 2, and 2 on the plasmid). These proteins are related 
to transposases found in several different families of inser- 
tion elements and in Tn3 family transposons. Thus, like most 
microbial genomes, the genome of 5. chlorophenolicum dis- 
plays evidence of continual onslaught by mobile genetic 
elements. However, the PCP degradation genes show no 
association with any of these elements. 

New genes that can help microbes survive and grow in 
the face of selective pressure from environmental toxins 
can also arise by gene duplication and divergence (Hughes 
1994; Bergthorsson et al. 2007). We carried out a Blast 
search of the 5. chlorophenolicum genome against itself 
to identify genes that are nearly identical and may have 
arisen by recent gene duplication. Only 38 proteins have 
>90% sequence identity to another protein in 5. chlorophe- 
nolicum. Of these, 13 are related to transposases; three 



of these genes are found in three identical copies scattered 
throughout the genome. Although the presence of highly 
similar transposase genes is not unexpected, there is little 
rhyme or reason to the identities and locations of the 
remaining duplicated genes. Seven genes annotated as 
encoding 2-hydroxychromene-2-carboxylate isomerase, 
a short-chain alcohol dehydrogenase, a hypothetical pro- 
tein, a nucleoside-diphosphate-sugar epimerase, an arabi- 
nose efflux permease, a Zn-dependent dipeptidase, and 
an outer membrane receptor protein are present between 
positions 667462 and 675176 on the — strand of chromo- 
some 2. A second copy of five of these genes are found 
in a cluster on the + strand of chromosome 2 between 
positions 1036969 and 1030812. Copies of the first two 
genes in the cluster are found together in a different loca- 
tion on chromosome 2. Nearly identical copies of the genes 
encoding the alpha and beta subunits of the E1 component 
of pyruvate dehydrogenase are found on chromosomes 1 
and 2 and nearly identical copies of a gene annotated as 
an acyl-CoA synthetase/AMP-acid ligase II are found in chro- 
mosome 2 but surrounded by completely different genes. 
The utility of these duplicated genes and the processes lead- 
ing to their location in distant parts of the genome are not 
clear. One possibility is that they were not actually dupli- 
cated in 5. chlorophenolicum but rather acquired by integra- 
tion into the genome of DNA fragments taken up from lysed 
cells of 5. chlorophenolicum or closely related Sphingomo- 
nads. Notably, none of the PCP degradation genes is found 
among the set of highly similar genes. 

Comparison of the Genomes of 5. chlorophenolicum 
and 5. japonicum 

The genomes of S. chlorophenolicum and S. japonicum are 
of similar size (4.57 and 4.46 Mbp, respectively) and contain 
a similar number of open reading frames (ORFs) (4, 1 59 and 
4,460, respectively). Both 5. chlorophenolicum and S. japo- 
nicum (Nagata et al. 2010) have a primary chromosome 
containing most of the genes for core processes (including 
glycolysis, the TCA cycle, amino acid and nucleotide biosyn- 
thesis, fatty acid oxidation, DNA replication, transcription, 
and translation) and a secondary chromosome. Sphin- 
gobium chlorophenolicum has a single plasmid (pSphChOI ), 
whereas S. japonicum has three (pUT1, pUT2, and pCHQ1). 

The third circle in each chromosome map in figure 2 
shows in green the positions of genes that encode proteins 
in 5. chlorophenolicum that have close homologs in 5. japo- 
nicum that exhibit >80% identity over >90% of the length 
of the 5. chlorophenolicum sequence. A total of 2,324 5. 
chlorophenolicum genes (1,931 on chromosome 1, 285 
on chromosome 2, and 108 on pSphChOI) have close ho- 
mologs in 5. japonicum (see supplementary table 2, Supple- 
mentary Material Online). In both chromosomes, regions 
with close homologs in 5. japonicum show a typical GC 
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Fig. 4. — Relationships between proteins encoded by (a) chromo- 
some 1, Ob) chromosome 2, and (c) the plasmid of Sphingobium 
chlorophenolicum and the best hits in the Sphingobium japonicum 
genome found by a Blast search using each 5. chlorophenolicum protein 
as a query sequence. The number of pairs is plotted as a function of % 
identity and % coverage of the query sequence for each replicon. (Only 
pairs for which E < 0.0001 are shown.) 



content of about 64%. Regions without close homologs 
(Sc islands) might have resulted from either loss of genes 
in 5. japonicum or acquisition of genes by HGT in 5. chlor- 
ophenolicum. Some of the Sc islands show a lower GC con- 
tent, suggesting that these regions may have been acquired 
by HGT. Most of the genes in the Sc islands are hypothetical 
proteins, but there are several predicted glycosy transferases 
(some predicted to be involved in cell wall biosynthesis), as 
well as some O-antigen ligases, some ABC transporters, 
cellobiose phosphorylase, a K+-transporting ATPase, and 
Type IV secretory pathway components. Notably, almost 
all the sequences associated with transposons and phage 
genes are found in Sc islands. 

A more detailed analysis of the relationships between 
proteins found in both 5. chlorophenolicum and 5. japoni- 
cum is shown in figure 4, which shows plots of sequence 
identity versus coverage for the top hit in the 5. japonicum 
genome for each 5. chlorophenolicum protein. (Coverage is 
defined as the length of the 5. chlorophenolicum query se- 
quence that is aligned to a sequence in 5. japonicum divided 
by the total length of the query sequence. The data were 
filtered to remove pairs for which the e value was 
>0.0001.) On this plot, homologs cluster in three regions: 
1) close homologs that share high sequence identity over 
most of the query sequence; 2) more distant homologs that 
share moderate sequence identity over most of the query 
sequence; and 3) homologs that share sequence identity on- 
ly over part of the query sequence. Notably, chromosome 1 
is highly enriched in close homologs, and chromosome 2 is 
modestly enriched in distant homologs. Together with 
the observation that most of the genes for core metabolic 
processes are present on chromosome 1, this observation 
suggests that chromosome 2 may preferentially collect 
horizontally transferred genes. 

Figure 4a shows that most of the proteins encoded on 
chromosome 1 share >80% identity with proteins in 5. ja- 
ponicum; indeed, 50% share >90% identity. We posit that 
close homologs with >80% identity were present in the 
most recent common ancestor of these two species; most 
are likely to be orthologs. The more distant homologs in re- 
gion 2 are unlikely to be orthologs derived from the most 
recent common ancestor, since they are much more diver- 
gent than the large number of close homologs in region 1 . 
The genes encoding these proteins may have been acquired 
by HGT independently from different sources in the two spe- 
cies; they may serve the same or different functions. Alter- 
natively, the 5. chlorophenolicum protein may indeed have 
derived from the most recent common ancestor, but the or- 
tholog in 5. japonicum may have been lost so that the best 
hit in the 5. japonicum genome is actually a paralog. Finally, 
there are 1 1 2 proteins for which significant sequence iden- 
tity is seen over only part of the query sequence (<75%); 
these include proteins in which conserved domains have 
been utilized in different structural contexts. 
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Fig. 5. — Scatter plot of gene conservation in chromosomes 1 and 2 of Sphingobium chlorophenolicum and Sphingobium japonicum. Red dots 
correspond to pairs of ORFs found when the two chromosomes are compared directly and blue dots to pairs of ORFs found when one genome is 
aligned to the reverse complement of the other. 



Although the dominant chromosomes of 5. chloropheno- 
licum and 5. japonicum share a common core of genes, 
there has been considerable rearrangement of genes since 
the common ancestor of these two bacteria. Figure 5 shows 
a scatter plot of gene conservation between the chromo- 
somes of 5. chlorophenolicum and 5. japonicum made using 
nucmer in the MUMmer package with the default param- 
eters (Delcher et al. 1 999, 2002; Kurtz et al. 2004). This plot 
shows a pattern known as an "x-alignment" (Eisen et al. 
2000) that is often seen in closely related bacteria. Many 
genes in chromosome 1 of 5. chlorophenolicum are found 



in comparable positions in chromosome 1 of 5. japonicum, 
as indicated by the red dots along the diagonal. However, 
a number are found in an inverted orientation (see blue 
dots) in positions that lie close to a diagonal perpendicular 
to the red diagonal. This pattern is believed to result from 
multiple inversions centered around either the origin or 
terminus of replication. There has evidently been consider- 
able remodeling of the genome via movement of blocks of 
genes within chromosome 1 since 5. chlorophenolicum and 
5. japonicum diverged from a common ancestor. The 
x-alignment pattern is less distinct for chromosome 2, 
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suggesting greater plasticity in this replicon. This is consis- 
tent with the lower number of genes for core processes 
on chromosome 2 described above. Furthermore, only 12 
of the 30 genes encoding core processes found on chromo- 
some 2 of 5. chlorophenolicum are found on chromosome 2 
of 5. japonicum, suggesting that movement of genes be- 
tween the two chromosomes has continued to occur since 
the divergence of 5. chlorophenolicum and 5. japonicum. 

Figure 6 depicts the correspondence between the posi- 
tions of very close homologs (>90% identity over >80% 
of the query length) in the entire genomes of 5. chlorophe- 
nolicum and S. japonicum. Light blue, orange, and dark blue 
lines connect the positions of genes in chromosome 1 , chro- 
mosome 2, and the plasmid, respectively, of 5. chlorophe- 
nolicum with the positions of homologs in the two 
chromosomes and three plasmids of 5. japonicum. As noted 
above, chromosome 1 in both species is densely populated 
with shared genes. A number of genes found on chromo- 
some 2 of 5. chlorophenolicum have homologs on chromo- 
some 1 of 5. japonicum, although the converse is not true. 
Notably, genes found in 5. chlorophenolicum but not in 5. 
japonicum are more heavily represented on chromosome 2 
than on Chromosome 1. Since chromosome 2 appears to 
carry many of the genes for degradation of organic com- 
pounds, this difference may be due to the availability of dif- 
ferent carbon sources in the environmental niches occupied 
by the two bacteria. Notably, the 5. chlorophenolicum plas- 
mid, pSPHCHOI, is comprised of a large region that is syn- 
tenic with a region of S. japonicum chromosome 1 (with the 
exception of a few small indels) and a smaller region that is 
syntenic with a region around the origin of 5. japonicum 
pCHQ1. This region encodes several proteins, including 
two chromosome partitioning proteins (ParA and ParB 
homologs) and the plasmid replication initiation protein 
(RepA). Thus, the 5. chlorophenolicum pSphChOI and 
S. japonicum pCHQ1 share an origin of replication and 
the associated genes, but the genes carried on the two plas- 
mids are not closely related. The two smallest plasmids in 
5. japonicum (pUT1 and pUT2) carry genes with no homo- 
logs in 5. chlorophenolicum. 

One hundred and eight of 125 genes on pSPHCHOI have 
>90% identity to genes on plasmid pSLGP in Sphingobium 
SYK-6 (a more distantly related Sphingomonad that is the 
closest relative of 5. chlorophenolicum and 5. japonicum 
for which a whole genome sequence is available), suggest- 
ing that these genes may have been present on a plasmid in 
the ancestor of 5. chlorophenolicum and 5. japonicum 
and may have been incorporated into chromosome 1 of 
5. japonicum after divergence of the two species. Move- 
ment of blocks of genes among the plasmids and chromo- 
somes in these organisms may have been facilitated by 
transposases, as TN3-family transposase elements are pres- 
ent in both plasmids and in the 5. japonicum chromosome 
near one end of the integrated region. 



Insights into the Origin of the PCP Degradation Genes 

The PCP degradation genes are found on chromosome 2. 
pcpA, pcpC, and pcpE are found in proximity to each other, 
although not within an operon, whereas pcpB and pcpD 
are found in an operon on the opposite side of the 
chromosome (see figs. 2 and 7). Expression of all the genes 
except pcpC (Orser et al. 1993) is induced in the presence 
of PCP, although it is not known whether the genes are 
regulated by PCP itself or by downstream products. The 
constitutive expression of pcpC suggests that a mechanism 
for its regulation by PCP has not yet arisen. 

In theory, the genes encoding the PCP degradation en- 
zymes might have originated by 1) recruitment of an ances- 
tral enzyme without gene duplication, requiring sharing of 
the enzyme between the new and original functions; 2) re- 
cruitment of a pre-existing enzyme, followed by duplication 
and divergence of the original gene to provide one copy that 
is specialized, although possibly not optimized, for its func- 
tion in PCP degradation; and 3) HGTof genes encoding ei- 
ther enzymes that have already specialized for their 
functions in PCP degradation or enzymes with inefficient 
promiscuous activities that became useful when 5. chloro- 
phenolicum encountered PCP. Comparison of the genomes 
of 5. chlorophenolicum and 5. japonicum allows us to 
determine which of these possibilities is most likely for each 
of the PCP degradation genes. If an enzyme originated by 
recruitment of an enzyme present in the most recent com- 
mon ancestor of S. chlorophenolicum and 5. japonicum, 
then a close homolog is likely to be found in 5. japonicum. 
If gene duplication and divergence has resulted in an 
enzyme specialized for its role in PCP degradation, we 
should find a close homolog in S. chlorophenolicum. If 
a gene has been acquired by HGT, we would likely find 
no close homologs in either 5. chlorophenolicum or S. japo- 
nicum and might find other signatures of HGT, as well. 

PCP hydroxylase (PcpB) and PcpD (TCBQ reductase) are 
encoded by an operon on chromosome 2 that also includes 
pcpR, which encodes a LysR-transcriptional regulator that 
has been reported to be essential for induction of expression 
of pcpB, pcpA, and pcpE (Cai and Xun 2002). PCP hydrox- 
ylase belongs to a family of flavin monooxygenases that hy- 
droxylate phenols. Such enzymes are common in soil 
bacteria because phenolic compounds derived from lignin 
provide an important source of carbon. The lack of a close 
homolog in 5. japonicum (see table 2) suggests that pcpB 
was not present in the most recent common ancestor of 
5. chlorophenolicum and 5. japonicum. Close homologs 
of pcpB are also not found in 5. wittichii, so it is unlikely that 
pcpB arose from an ancestral gene that was lost in 5. japo- 
nicum. The low pairwise sequence identities between PCP 
hydroxylase and other flavin monooxygenases in 5. chloro- 
phenolicum (see table 2) suggest that pcpB did not arise by 
recent duplication and divergence of a previously existing 
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Fig. 6. — Plot of gene mappings between replicons of Sphingobium chlorophenolicum (blue, orange, and dark blue) and Sphingobium japonicum 
(green). Blue, orange, and dark blue lines connect the positions of genes in chromosome 1, chromosome 2, and the plasmid, respectively, of 5. 
chlorophenolicum with the positions of homologs in 5. japonicum with >80% identity over >90% of the S. chlorophenolicum sequence. Diagram was 
generated using Circos (Krzywinski et al. 2009). 



5. chlorophenolicum gene. These findings, along with the 
observation that pcpB and pcpD fall in a region of chromo- 
some 2 with relatively low GC content (see fig. 2), suggest 
that pcpB was acquired by HGT from an unknown source. 
This hypothesis is supported by the observation that en- 
zymes with 72-98% identity to PcpB are found in Novos- 



phingobium lentum (Tiirola, Mannisto, et al. 2002), 
several polychlorophenol-degrading Sphingomonads from 
Finland (Tiirola, Wang, et al. 2002), Sphingomonas sp. 
UG30 (Cassidy et al. 1999), and in several uncultured bac- 
teria from environmental samples collected from PCP- 
contaminated soils (Beaulieu et al. 2000). The high sequence 
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Fig. 7. — Gene neighborhoods surrounding the PCP degradation 
genes. pcpR and pcpM encode transcriptional regulators that control 
expression of pcpBD, pcpA, and pcpE. (a) 7, predicted transcriptional 
regulator; 2 putative NADP-dependent oxidoreductase; 3, hypothetical 
protein; 4, glycosyl transferase family 2; 5, outer membrane receptor 
protein, mostly Fe transport; 6, outer membrane receptor for ferrienter- 
ochelin and colicins; 7, glucose/sorbosone dehydrogenase; (b) 8, outer 
membrane cobalamin acceptor protein; 9, predicted esterase; 10, 
formyltetrahydrofolate deformylase; 11, methenyltetrahydrofolate cyclo- 
hydrolase/methylene tetrahydrofolate dehydrogenase; 12, Kef-type K + 
transport systems, membrane components; 13, metal-dependent hydro- 
lase, (3-lactamase superfamily III; 14, outer membrane receptor proteins, 
mostly Fe transport; 15, predicted glutathione 5-transferase; 16, putative 
LysR-type transcriptional regulator. 



identities between these genes suggest that pcpB has 
been transferred among a number of bacteria in PCP- 
contaminated environments, a conclusion previously reached 
by Tiirola, Wang, et al. (2002). Strikingly, PcpB has <35% 
identity to most flavin monooxygenases from non-PCP 
degrading bacteria. The high level of sequence divergence 
makes it difficult to discern what the original substrate for 
PCP hydroxylase might have been before its recruitment into 



the PCP degradation pathway. (Careless annotation transfers 
have led to a proliferation of genes that are annotated as 
encoding PCP hydroxylase. The proteins encoded by these 
genes are homologous to PcpB, but have only about 35% 
identity to PcpB, a level at which substrate specificity cannot 
be confidently assigned. Furthermore, the ability of the cor- 
responding microbes to degrade PCP has not been assessed. 
These assignments are probably wrong in the majority of 
cases.) 

TCBQ reductase (PcpD) catalyzes the second step in the 
PCP degradation pathway (Dai et al. 2003). Reduction of 
a benzoquinone is an uncommon step in aromatic degrada- 
tion pathways. Benzoquinones are formed as a result of hy- 
droxylation of a phenol only when there is a leaving group, 
such as a chlorine or nitro group, at the position of hydrox- 
ylation. Most commonly, there is a hydrogen at the position 
of hydroxylation, and the product is a hydroquinone (see 
fig. 8). 

TCBQ reductase is not closely related to any other enzyme 
in 5. chlorophenolicum or 5. japonicum (see table 2), 
suggesting that it, like pcpB, was acquired by HGT. Enzymes 
that carry out comparable reductions of benzoquinones 
produced by phenol monooxygenases are found in three 
other microbes, Pseudomonas sp. Strain WBC-3 (Zhang 
et al. 2009), Alcaligenes sp. strain NyZ215 (Xiao et al. 
2007), and Cupriavidus necator J MP 1 34 (formerly Ralstonia 
eutropha) (Belchik and Xun 2008), which degrade p- 
nitrophenol, o-nitrophenol, and trichlorophenol, respec- 
tively. Notably, these three benzoquinone reductases and 
TCBQ reductase belong to three different superfamilies 
(see fig. 9). The recruitment of enzymes from three different 
superfamilies to catalyze reduction of quinones in pathways 
for degradation of four different anthropogenic compounds 
suggests that, in each case, a bacterium has taken 



Table 2 

Closest Homologs of PCP Degradation Genes in Sphingobium chlorophenolicum, Sphingobium japonicum, and Sphingobium sp. SYK-6 



Closest Homolog in Closest Homolog in Sphingobium sp. 

5. chlorophenolicum Closest Homolog in 5. japonicum SYK-6 





Protein 


% Identity 


Protein 


% Identity 


Protein 


% Identity 


PcpB (gi2 1362858) 


FAD-binding 
monooxygenase 
(gi 10740769) 


28 


FAD-binding 
monooxygenase 
(gi2941 46540) 


29 


p-hydroxy-benzoate 

hydroxylase (gi347528071) 


21 


PcpD (gi2 1362854) 


Phthalate 4,5-dioxygenase 
subunit (gi 1 0723764) 


42 


Vanillate monooxygenase 
oxidoreductase subunit 
(gi294023794) 


42 


Putative flavodoxin 
reductase (gi347430572) 


27 


PcpC (gi 22417110) 


Glutathione S-transferase 
domain-containing 
protein (gi 10724309) 


25 


2,5-dichlorohydro-quinone 
dehalogenase 
(LinD) (gi294023795) 


25 


Glutathione S-transferase 
(gi347527108) 


25 


PcpA (gi 3760223) 


None 




2,6-dichlorohydro-quinone 
1,2-dioxygenase (LinEb) 
(gi294146914) 


93 


None 




PcpE (gi 22417101) 


Maleylacetate 
reductase (gi334342570) 


53 


Maleylacetate reductase (LinF) 
(gi2941 46911) 


91 


Alcohol dehydrogenase 
(gi347529189) 


29 
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Fig. 8. — Products formed by hydroxylation at a position bearing a hydrogen (a) and a chlorine (b). 



advantage of a different pre-existing protein when faced 
with the need to degrade an aromatic compound that car- 
ries a leaving group at the site of the initial hydroxylation 
reaction. 



While PC P hydroxylase catalyzes a reaction similar to that 
of its homologs, the progenitor of TCBQ reductase may have 
served a different function. Although there are several fam- 
ilies of known quinone reductases, TCBQ reductase is not 
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Fig. 9. — Benzoquinone reductases have been recruited from three different superfamilies in microorganisms degrading PCP {Sphingobium 
chlorophenolicum), 2,4,6-trichlorophenol {Cupriavidus necator JMP134), p-nitrophenol {Pseudomonas sp. Strain WBC-3), and o-nitrophenol 
(Alcaligenes sp. strain NyZ21 5). The benzoquinones shown (from top to bottom) result from hydroxylation of a phenol at a position bearing a chlorine 
atom (for PCP and 2,4,6-trichlorophenol) or a nitro group (for p- and o-nitrophenol). 
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related to any of these. Rather, TCBQ reductase is most 
closely related to proteins that serve as the reductase com- 
ponent of two-component dioxygenases that initiate aero- 
bic degradation of aromatic compounds that lack hydroxyl 
or amine substituents (Dai et al. 2003). These proteins con- 
tain a flavin that accepts electrons from NAD(P)H and trans- 
fers them one at a time to an iron-sulfur cluster, which then 
donates electrons to the active site of the oxygenase com- 
ponent of the enzyme. TCBQ reductase is also more distantly 
related to ferredoxins that restore activity of extradiol dioxy- 
genases that have been inactivated during substrate turn- 
over (Polissi and Harayama 1993; Hugo et al. 1998, 2000; 
Tropel et al. 2002). (This inactivation occurs after 0 2 reacts 
with Fe(ll) at the active site of the enzyme-substrate complex 
to form superoxide and Fe(lll). If the superoxide diffuses out 
of the active site before reacting with the substrate, the iron 
atom is left in the inactive Fe(lll) state.) The reaction catalyzed 
by TCBQ reductase could certainly take advantage of this 
protein architecture; donation of hydride to a flavin, fol- 
lowed by transfer of single electrons to the iron-sulfur cluster 
and thence to bound TCBQ would form first a semiquinone 
and then TCHQ after transfer of a second electron. Thus, ad- 
aptation of the progenitor protein to this new role may have 
required little more than enhancement of a small molecule- 
binding site adjacent to the iron-sulfur cluster. 

TCHQ dehalogenase (PcpC) catalyzes the third and 
fourth steps in the PCP degradation pathway. Each reductive 
dehalogenation step results in oxidation of two molecules 
of glutathione to glutathione disulfide. Similar reductive 
dehalogenation steps occur during lindane degradation in 
5. japonicum (Miyauchi et al. 1 998) and during degradation 
of PCBs in Burkholderia xenovorans LB400 (Batels et al. 
1999; Tocheva et al. 2006). All of these dehalogenases 
are members of the glutathione 5-transferase family. Most 
enzymes in this superfamily catalyze the nucleophilic attack 
of glutathione upon an electrophilic substrate to form a glu- 
tathione conjugate. However, a few, such as TCHQ dehalo- 
genase (Warner et al. 2005; Warner and Copley 2007b), 
maleylacetoacetate isomerase (Polekhina et al. 2001), and 
maleylpyruvate isomerase (Marsh et al. 2008), catalyze 
more complicated reactions in which additional steps occur 
before and/or after the canonical attack of glutathione upon 
an electrophilic intermediate. 

Although there are numerous members of the glutathi- 
one S-transferase superfamily in 5. chlorophenolicum, 
TCHQ dehalogenase has no more than 25% identity to 
any of them, so pcpC likely did not arise by divergence from 
another gene in 5. chlorophenolicum. Close homologs 
of pcpC are not found in the genome of 5. japonicum, so 
pcpC was probably not present in the ancestor of 
5. chlorophenolicum and 5. japonicum. The GC content 
of pcpC is 61 %, and it falls in a region of chromosome 2 
with unusually low GC content. Thus, pcpC was most likely 
acquired by HGT Strikingly, 5. japonicum LinD, which cata- 



lyzes the reductive dehalogenation of 2,5-dichlorohydroqui- 
none to chlorohydroquinone (Miyauchi et al. 1 998), has only 
25% identity to TCHQ dehalogenase. TCHQ dehalogenase 
and LinD are the only known reductive dehalogenases in 
aerobic bacteria that act on aromatic compounds. The high 
divergence between these enzymes suggests that they were 
recruited independently from different members of the GST 
superfamily, possibly due to selective pressure to degrade 
PCP and lindane, respectively. The function of the progen- 
itors of these enzymes is not clear. Residues in the active sites 
of TCHQ dehalogenase and LinD resemble those in the ac- 
tive sites of some maleypyruvate isomerases, and TCHQ de- 
halogenase has low activity with maleylacetone, an 
analogue of maleylpyruvate (Anandarajah et al. 2000). Pro- 
miscuous activities can provide clues to the activity of an an- 
cestral enzyme, so TCHQ dehalogenase may have originated 
from a maleylpyruvate isomerase. However, this cannot be 
stated with certainty. 

2,6-Dichlorohydroquinone dioxygenase (PcpA) catalyzes 
the cleavage of 2,6-dichlorohydroquinone (DCHQ). DCHQ 
dioxygenase is 93% identical to LinEb in 5. japonicum, 
which also cleaves DCHQ. (Despite its name, LinEb is not 
involved in degradation of lindane [Endo et al. 2005].) Thus, 
this enzyme was likely present in the most recent common 
ancestor of 5. chlorophenolicum and S. japonicum. Whether 
cleavage of DCHQ is the natural function of this enzyme is 
not clear, as studies of its substrate range and catalytic 
efficiency have not been reported beyond an early report 
that shows that hydroquinone and chlorohydroquinone also 
serve as substrates (Ohtsubo et al. 1999). Sphingobium 
japonicum does not degrade PCP, so the role of LinEb in this 
bacterium is not known. {Sphingobium japonicum contains 
a homolog of LinEb called LinE that cleaves chlorohydroqui- 
none in the lindane degradation pathway. LinE and LinEb 
share 53% identity and DCHQ dioxygenase and LinE share 
52% identity.) 

The final two steps in degradation of PCP are catalyzed by 
maleylacetate reductase (PcpE), which carries out the succes- 
sive reductions of 2-chloromaleylacetate to maleylacetate 
and of maleylacetate to (3-ketoadipate. PcpE has 91% 
identity to maleylacetate reductase from 5. japonicum (LinF), 
indicating that maleylacetate reductase was found in the 
progenitor of 5. chlorophenolicum and 5. japonicum. Thus, 
PcpE has been probably been recruited to carry out the reduc- 
tion of 2-chloromaleylacetate in addition to the more typical 
reduction of maleylacetate. The 5. chlorophenolicum ge- 
nome also encodes two other maleylacetate reductase 
homologs, Sc_00030750 and Sc_00027830, which have 
53% and 55% pairwise sequence identities to PcpE, respec- 
tively. These enzymes have not been characterized, but one 
or both may play a role in reduction of 2-chloromaleylacetate 
based upon the observation that substantial 2-chloromaley- 
lacetate reductase activity is present in a strain in which pcpE 
has been knocked out (Cai and Xun 2002). 
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Conclusions 

Comparison of the genome sequences of 5. chlorophenoli- 
cum and 5. japonicum suggests that the first three genes in 
the PCP degradation pathway, pcpB, pcpD, and pcpG, were 
acquired by HGT. Since these genes were incorporated into 
two different places in the genome, at least two different 
HGT were involved. The documented spread of nearly iden- 
tical pcpB genes among bacteria in PCP-contaminated 
soils suggests that transfer of pcpB, and possibly pcpD as 
well, may have occurred within the last century since the 
introduction of PCP into the environments. 

The last two genes in the PCP degradation pathway, pcpA 
and pcpE, were inherited from the progenitor of 5. chloro- 
phenolicum and 5. japonicum. Pathways for degradation of 
aromatic compounds typically begin with specialized 
enzymes that convert the initial compounds to catechol 
or hydroquinone intermediates that are then cleaved by 
intradiol or extradiol dioxygenases. (Intradiol dioxygenases 
cleave catechols between the two hydroxyl groups. Extradiol 
dioxygenases cleave catechols or hydroquinones adjacent to 
one of the hydroxyl groups.) The ring cleavage products are 
then processed to C3 and C4 carboxylates that feed into 
central metabolism. Thus, numerous "upper" pathways fun- 
nel into a small number of "lower" pathways. The enzymes 
encoded by pcpA and pcpE likely existed as part of a stan- 
dard lower pathway for degradation of naturally occurring 
compounds. The close proximity of pcpA and pcpE in the 
genome is consistent with this interpretation. 

Duplication and divergence is believed to be the most 
common mechanism for evolution of enzymes with novel 
capabilities. This process allows exploration of sequence 
space while the original activity of the encoded enzyme is 
provided by a copy that deviates little, if at all, from the orig- 
inal sequence. Ultimately, extraneous copies of the gene are 
eliminated, leaving one copy encoding the original enzyme 
and one copy encoding the novel enzyme. Analysis of the 
5. chlorophenolicum genome sequence does not reveal 
evidence for recent duplication and divergence of any of 
the PCP degradation genes within S. chlorophenolicum itself. 
Duplication and divergence may have occurred in the un- 
known soil bacterium in which the genes originated, and only 
the genes encoding enzymes useful for degradation of PCP 
may have been transferred to S. chlorophenolicum and other 
PCP-degrading Sphingomonads. 

Gene amplification is commonly observed when ineffi- 
cient enzymatic activity limits the ability of a bacterium to 
degrade a carbon source. For example, the benzoate deg- 
radation genes in a strain of Acinetobacter AD P-1 that lacks 
the transcriptional regulators that normally induce these 
genes are amplified up to 20-fold when benzoate was sup- 
plied as a sole carbon source (Reams and Neidle 2003). 
Given the inefficiency of the initial enzymes in the PCP deg- 
radation pathway, it is intriguing that there is no evidence 



for amplification of genomic segments containing the genes 
encoding these enzymes. In this case, cells carrying dupli- 
cated or amplified genes may be at a disadvantage. Ampli- 
fication of a segment of the genome carrying pcpB and 
pcpD would indeed increase flux through the initial two 
steps but would probably not improve fitness. The waste 
of NADPH and the concomitant production of H 2 0 2 by 
the uncoupled PCP hydroxylase reaction would increase 
in proportion to the enzyme concentration. Additionally, 
the increased concentration of TCHQ formed would inhibit 
TCHQ dehalogenase, leading to a buildup of TCHQ, which 
would be expected to cause redox cycling and depletion of 
cellular reductants. On the other hand, duplication of a seg- 
ment of the genome carrying pcpC would provide no selec- 
tive advantage because the initial step, conversion of PCP to 
TCBQ, limits the flux through the pathway (McCarthy et al. 
1997). Thus, duplication of this region might simply cause 
a metabolic burden without providing any benefit. If pcpB, 
pcpD, and pcpC were colocalized, then duplication of all the 
genes together might be beneficial. However, given the in- 
corporation of pcpBD and pcpC into two distant regions of 
chromosome 2, gene duplication does not provide a path- 
way for further evolution at this point. 

Supplementary Material 

Supplementary tables 1 and 2 are available at Genome Bio- 
logy and Evolution online (http://www.gbe. oxfordjournals. 
org/). 
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